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ON ■ Abstract 
^ . 

For the time development of a single system in the quantum jump approach 
O I or for quantum trajectories one requires the conditional (reduced) Hamiltonian be- 

tween jumps and the reset operator after a jump. Explicit expressions for them are 
derived for a general A^-level system by employing the same assumptions as in the 
usual derivation of the Bloch equations. We discuss a possible minor problem with 
positivity for these expressions as well as for the corresponding Bloch equations. 
PACS numbers: 42.50, 32.90 +a 
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1. Introduction 



In order to better understand photon statistics and photon counting processes for single 
atoms, in particular the macroscopic dark periods ("electron shelving") in the Dehmelt 
V system the quantum jump approach was developed in Gottingen prior to 1991 
S ■ 01 5 consequently analytic aspects were stressed. The independently proposed 

Monte-Carlo wave function approach is essentially equivalent, but stressed more the 
computational aspects. The approach by quantum trajectories H] was also independently 



■ proposed and is also essentially equivalent. There are various papers which can be re- 
^ ■ garded as precursors to the quantum jump approach, in particular Refs. [§, ^. Some 



related work P, Rl, |10| on photon statistics will be discussed in Section 5. 



Our analytic methods have been applied to interference effects in single atoms [IT], |12 



T^. Recently, the quantum jump approach has been extended to study conditional 



fluorescence spectra, such as the spectrum of an atom in an extended light period flS 



In the derivation of the quantum jump approach one envisages photon measurements 
in rapid succession at times At apart. If a photon is detected at some time it is, for 
simplicity, assumed to be absorbed. This assumption, however, is not essential [|l6l. At 



should be much smaller than the level life-times, but it cannot be too small because then 



one would encounter the quantum Zeno effect jl^ ; therefore At should be larger than the 
inverse optical frequencies. A reasonable range is At ^ 10~^^s. 

Between broadband photon detections the time development of an atom is then de- 
scribed by a reduced, or rather conditional, Hamiltonian -ffcond, the condition being that 
no photons are found. It is non-hermitian and allows the calculation of the detection 
or emission probability. Once a photon is detected, the atomic state has to be reset (a 
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"jump"), usually to the ground state, but for general A^-level systems the reset states 
may be quite complicated . Thus the state of a single atom is given by a random path 
("quantum trajectory") determined by the conditional Hamiltonian -ffcond and the reset 
states. 

In Ref. one of us derived a simple expression for the reset state in the two limiting 
cases of either widely separated or closely neighbored transition frequencies. The general 
case seemed to require much smaller At's, which are difficult to justify physically. 

The main purpose of this paper is to re-investigate the reset state using the same 
methods and approximations as in the derivations of the Bloch equations |T^, |^, mainly 
the Markov approximation. It will be shown that in the general case the (super-)operator 
which resets the state may depend in a complicated oscillatory way on the choice of 
At. However, this "physical" reset operator can be replaced by a simpler refined or 
idealized reset operator which gives the same results for all photon counting processes. 
This idealized reset operator turns out to have the same form as the one obtained in Ref. 
by using very small At's. This is explained in detail in Sections 3, while in Section 
2 a short derivation of the conditional (or reduced) Hamilton is given. In Section 4 it 
is shown that the collection of quantum trajectories satisfies the Bloch equations, both 
for the physical and for the idealized reset operator. In the last section we discuss our 
results and the positivity of the expressions. It turns out that in some cases the reset 
operator and the general damping matrix in the conditional Hamiltonian can have a small 
negative part. A related problem also occurs in Bloch equations where in some examples 
the solution develops a small negative part for very small times and thus temporarily 
looses its positivity But this is probably of no practical relevance. 



2. The conditional Hamiltonian Hcond for a general N-level system 



With an external field Ee(t), the Hamiltonian for complete system, atom plus quantized 
electromagnetic field E = E^"*") + is given in the Schrodinger picture, in the limit of 
long wave-lengths and in the rotating-wave approximation, by |20| 



H = Y.^^i K) (^1 +Hf+ [eD(-) • Ei+\t) + h.c] + [eD^") ■ E^ + h.c. 
= + Hp + HAL(t) + Haf 



where 



D(-) 
E(+) 



i>j 

kA 



1/2 



2eoV 



(2) 



and a frequency cutoff is included. V is the quantization volume, later taken in infinity, 
and i > j means Ui > Uj. 

Now one imagines that photon measurements are performed at times ti = At, ^2 = 
2 At, ■ ■ ■ on an ensemble of systems which, at time t = 0, is supposed to be in the initial 
state \Oph)\'ip)- By the von Neumann - Liiders projection postulate ||2T[ the subensemble 
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for which no photons are detected until time t„ is described, with Pq = \0ph)lA{0ph\, by 

PoUitn,t.-At)Po---PoUih,0)\OpM) = Kh)U,onditn,Om • (3) 

The norm squared of this is the probabihty of finding no photons for the measurements 
between and tn- t/^cond gives the time development of an atom under the condition that 
no photon is observed until time t, and it will now be determined explicitly by simple 
second order perturbation theory. Going over to the interaction picture with respect to 
H% + H^p one has 

H,{t)=Hi^{t) + H'^p{t) (4) 

which is obtained by replacing and akA in the original interaction Hamiltonian by 

|i)(j|e*'^*-'* and OkAe"*'^''*, respectively, with 



^ij = — ■ 



We now calculate, for ti <t' < tj+i, 

{%h\j^Ui{t\U)\%n) . (5) 

In the first-order contribution only Hj^]^{t) remains since 

{Oph\H'^Fit)Kh) = . (6) 

The second order is, by Eq. (H), 

- h-' f dt" {{%n\H^^p{t')Hip{t")\%n) + H^AL{t')H'AL{t")] . (7) 

Jti 

Now, if the external field Ee(t) is smooth in time, e.g. given by laser, and not wildly 
fluctuating like a thermal or chaotic field, then the second part in Eq. (|^) contributes 
a term of higher order in At and can therefore be omitted. For a chaotic external field, 
however, this part may give rise to a contribution of the same order and then has to be 
retained. In this case one can no longer work with state vectors ("wave functions") but 
has to use (conditional) density matrices. A particular example of this is treated in Ref. 



12|. 

Thus, supposing a smooth external field, the second-order contribution becomes 



n-' f dt" E K)0-H(^lE$^(D^.-^kA)(ekA-D 



ijlm kA 
i>j, l>m. 



•'0 kA ^^0»/ 

One can now use properties of the correlation function 
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K,,em{r) = E 7p!^i^J^ ■ ^kA)(ekA " ^im) e-^^^"^— . (9) 

kA "^^onV 
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With = A^k/{27i)^ one can perform the hmit V oo, and the sum over k becomes 
an integral over u, with a suitable frequency cutoff, and an integral over the unit sphere. 
The correlation function has an effective width of the order of around r = 0, and for 
t' — ti 3> uj^j^ one can therefore extend the r integration in Eq. (^) to infinity . This 
is equivalent to the approximation 

dre^(^*--^>')" ^ TT^K - cu£^) + (10) 

t^k — ^im 

and corresponds to the usual Markov approximation in the derivation of the Bloch equa- 
tions [p!8| , For the second-order contribution one then obtains 

- E K)(^|e*""*' / d^'k . \ D,, ■ J2 gkAgkA • D,,(7r(5(a;k - ooe^) +tV ) . (11) 

i.l>3 



The principal- value term is analogous to a level shift and is often omitted The 
last integral equals Tijj^ where 

^ijM = rrDij ■ Tikt\i^ktf + principal value term (12) 

Hence, integrating over t' from ti to tj+i and using 1 + 5 ~ e'' for small 5, we obtain 



{Oph\Ui{U+i,U)\Oph) = exp |-^^ dt' I^HMf) - E r.,,-^e*"«*'K)(^||| • (13) 
For small At this can be replaced by a time-ordered exponential and thus, with t = tn, 

11 MUiiu, u_,)\Oph) = r exp dt' HAL{t') - E r.,,-^e^"^'*'K)(^| (14) 



where the product sign on the l.h.s. includes an ordering in an obvious way. 
Since 

(Op^|f/(t„t,„i)|Op;.) = e-^^AW?i(Op^|f/,(t.,t,_,)|Op,)e^^-V-i/'* (15) 
and since, for t = tn = nAt, 

n 

f/cond(t,0) = l[{Oph\U{ti,ti_i)\Oph), 
1 

we obtain, on a coarse-grained time scale, from Eqs. (|T^ and (^) 

f/cond(t, 0) = r exp rft' + i/^i(t') - 2^1 X; r,,,,|z) (£| I I (16) 

which is the transformation of Eq. (|l^) back to the Schrodinger picture. 
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Thus, with the atomic operator T defined as 



r= Er.,,.K)(^l (17) 



it] 

i,e>j 



the conditional Hamiltonian for an iV-level atom with no photon emission until time t is, 
on the coarse-grained time scale, given by 

i^cond(t) =K + HMt) - . (18) 

For initial atomic state the probability to find no photon until time t is thus given 

by 

||f/cond(t,0)|7/')||2 , (19) 

and the probability to find the first photon in (t, t + At) is the difference of this expression 
for t and t + At. Thus, on the coarse grained time scale the probability density w(t) for 
the first photon is the negative derivative of Eq. (pIS]), 



<r) = -|||f/_,(t,0)|^)||2 

= (v^ir + riv^). (20) 

When one lets At become smaller and smaller the above derivation shows very nicely 
how and where the quantum Zeno effect [|l7l turns up in a very natural way.If At is chosen 
much smaller than the inverse optical frequencies, the last exponential in Eq. (|]) can be 
replaced by 1, and the integral becomes proportional to t' — tj.Eq. (|13D is then replaced 
by 

{Oph\UiiU+uti)\Oph) = exp I I HALiU)At - ih const ^ ri,-,,e''^«*'|z)(£|(At)2 

V I i,l>] 

(21) 

The product of these operators then becomes, for At — > 0, 



Texp|-^|Jrft'{^AL(t')}} • (22) 



This is a purely atomic operator, and hence the time development of the field becomes 
frozen, i. e. for At ^ one always remains in the vacuum. For this reason one cannot 
choose At arbitrarily small in the quantum jump approach. 



3. The reset operator 

In this section we determine the state (or density matrix) of an atom after a broadband 
detection of a photon, under the condition that shortly before no photon was found. We 
therefore consider an ensemble where no photon are present at time t„. The ensemble is 
thus described by p(t„) = \Qph) pA(tn){^ph\ where pA is the density matrix of the atoms. 
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The state of the subensemble for which photons are detected by a non-absorptive mea- 
surement at time tn+i = tn + is given, in view of the von Neumann-Liiders projection 
postulate ||21|, by 

Pip(tn+i)Pi/tr(-) (23) 

where 

Pl = l-|0p;,)l^(0p;,| . (24) 

Note that Eq. ( p3D still contains the photons. 

After a photon measurement by absorption no photons are present any longer and it 
was argued in Ref. that the resulting reset state is obtained from Eq. ( PB] ) by a partial 
trace over the photons, i.e. by 

Kh) {tl^hPlp{tn+l)Pl) {%h\/tT{-) . (25) 

The physical reason for this is that for the atomic description alone it should make no 
difference in infinite space whether or not the photons are absorbed, as long as they are 
sufficiently far away from the atom and no longer interacting with it [Q. Eq. ( [25| ) can 



be calculated by perturbation theory for [//(t„+i, t„), as in Section 2. Now, the first-order 
contribution suffices and one obtains in a straightforward way for the reset state 

7^(A^)pA(^n) = iTphPlU{tn+l,tn)\%h)pA{tn){%h\U{tn+l,tnyPl 

= e-'^'A^"^ E \j){i\pA{tnm{m\e^">'/'' 

ijim 



I dt' I dt"e->^(*'-*")— E 7^(D,, ■ ekA)(^kA ■ D,^) . 
Jo Jo ^ iSohV 

(26) 



Note that the atomic trace tr^7?.(At)p^(tji) gives the probability for a photon to be 
found at time tn+i under the condition that no photons were found at time t„. 

The external field drops out in the first-order contribution since its action during the 
short time At is of second order only. To apply the Markov property, we decompose the 
rectangular integration domain over t' and t" in Eq. (EB) into two triangles, leading to 



25| 



At ct' rAt ft" 



JO JO JO 



(27) 

As in Eqs. (||) and (p!0| ) the inner integrals can be replaced by 716(00^ — ujim) and 7r(5(a;k — 
uiij), respectively, plus principal values. In the limit ^ 00 the sum over k becomes an 
integral over d^k as in Eq. (11). With Tijj^g given by Eq. (12) we thus obtain 



ij£m 



(28) 

Up to normalization this is the state of an atom after a photon detection at time 
tn + At, under the condition that no photons were found at time Its trace gives the 
probability for this event. 
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Even after normalization the state will in general depend in an oscillatory way on At, 
except in two limiting cases which were discussed in Ref. 0|. If two optical transition 
frequencies Uij and uim are far apart, then \uij —u!£m\At 3> 1 and the integral in Eq. ( ]28|) 
vanishes. The other case is when two optical transition frequencies are very close so that 
\uJii —uJirn\At <^ 1. Then the integral in Eq. (^) is essentially 1, and in this case one may 



obtain interesting coherence effects which were discussed for the A system in Ref. |13 



There is a close connection between the reset operator, the conditional Hamiltonian, 
and the Bloch equations. By Eq. (|2^) one can write Eq. (E^) for the reset state as 



(29) 

with p(tn) = \Oph)pA(tn){Oph\- The first term on the right hand side is the definition 
of the atomic density matrix in the Bloch equations, while the second is given by the 
conditional time development operator. Hence one can obtain the reset operator also 
from a knowledge of the Bloch equations and -ffcond- 

For the general case the reset operator in Eq. ( ^8]) is cumbersome to work with, and 
it will now be shown that one can replace it by a refined - and simpler - expression which 
gives equivalent results for all photon counting questions. 

We define the atomic superoperator J by 



JpA= J2 {'^jiem + '^emji}\j){i\pA\^){fn\ . (30) 

ij£m 

Formally one has 

J = 7^(A^)/A^ (31) 

but of course this limit is not really physically allowed. However, the photon counting 
probabilities obtained from J for physically allowed At's turn out to be in agreement with 
those obtained from TZ{At). 

As in Ref. we define the atomic superoperator S{t,tQ) by 

S(t,tQ)pA = Ucond(t,'to)PAUcond(t,to)* ■ (32) 

Based on J as reset operator the probability density w{ti, ■ ■ ■ ,Tk; [0, t]) for finding a photon 
exactly at the times ri, ■ ■ ■ , in [0, t] is given by [|| 

w{ti, ■■■,Tk; [0,t]) = ti a{S (t, Tk)J S {Tk, Tk-i)J ■ ■ ■ JS'(ri,Oi)pA(0)) (33) 

To see how this compares with the photon counting probabilities obtained from TZ{At) 
we first express the latter by means of J as follows. Let To(t, to) be the superoperator for 
the free time evolution of atomic density matrices, 

and let T(t, tg) be the time evolution including driving and damping, i.e. given by the 
solution of the Bloch equations. With Tq Eq. (^) can obviously be written as 

/•At 

7^(A^) = / dt'To{At-t',0)JTo{t',0) . (35) 
Jo 
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We now use this expression for TZ{At) to calculate the probability of finding a photon 
in the small time interval [ti, ti+i] and none between and tj and none between tj+i and 
t. With a simple change of integration range in Eq. ( pS]) the desired probability is given 
by 

tTAS{t,U+i)']Z{At)S{U,0)pA{0) = r^'rfritr^^(t,t,+i)fo(t,+i,ri)jfo(ri,t,)5(t,,0)p^(0) . 

Jti 

(36) 

Now, since At = tn+i — tn is small compared to the inverse damping and driving, one may 

- during At - replace Tq by T or by S. Here we replace it by S and obtain for Eq. (p^) 

dntTA{S{t, t,)JS{t,, O)p^(O)) . (37) 

Jti 

The integrand is just the probability density w{ti; [0,t]) of Eq. ( p3D associated with J 
for k = 1. It is clear that this argument immediately carries over to more than a single 
photon detection. The essential point is seen to be the fact that during a single time 
interval At the damping has negligible effect. 

We thus have shown that for all photon countings we can use the refined - or idealized 

- reset operator J and the densities of Eq. (|33|). 



4. Connection with Bloch Equations 

In Ref. it had been shown that an ensemble of atoms whose individual time devel- 
opment proceeded according to -ffcond and the reset operator J after photon detections 
("jumps") also obeys the Bloch equations. If instead one uses the reset operator lZ{At) 
from Eq. (pSj), to what extent do the Bloch equations still hold? 

The subensemble of atoms with no photons until t is described by S'(t, O)pa(O) and 
the subensemble with last photon detection before t at time ti by 

S{t,ti)iZ{At)pA{t,^i) 

where t = nAt and ti = iAt. Thus the complete ensemble of atoms is described at time 
thy 

PAit) = sit,o)pAio) + J2 s{t,ti)n{At)pA{ti^i) . 

i=l 

With Eq. (|35|) and a change of integration variable one obtains 

PA{t) = S{t, O)pa(O) + E / drSit, t,)To(t„ t)JTo{t, t,_^)pA{t,-i) . (38) 

As noticed before, during a single time interval At the driving and damping has negligible 
effect, and therefore one can replace in Eq. (p8|), for this short time interval, To(ti, r) by 
S(ti,T) and To(r, tj_i) by T(r, tj.i), where the latter is the time development operator 
of the Bloch equations including driving and damping. Since T(r, tj_i)p^(tj_i) = Pa^t), 
one obtains with this approximation 

n— 1 /•t-f.i 

PA{t) = S{t, O)pa(O) + T. I drSit, T)JpA{r) 

i=l •'^i 
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and thus 

PA{t) = S{t, 0)pa{0) + f drSit, r)JpA{r) (39) 

which is just the equation one would have obtained with J as reset operator. Introducing 
a coarse-grained time scale one can now differentiate with respect to t and obtains 

PA = -J^ {HconApAit) - PA(t)^^cond} + J PAit) 

which agrees with the usual Bloch equations p. 



5. Discussion 



The idea of the quantum jump approach is to describe the radiating atom between 
photon detections by a conditional (or reduced) time evolution operator giving the time 
development under the condition that no photon has been detected. After a photon 
detection one has to reset the atom to the reset state ("jump"), with ensuing reduced 
time development, and so on. For a driven system with many emissions one then obtains 
a stochastic path ("quantum trajectory"). An ensemble of such paths for an ensemble of 
driven atoms satisfies the Bloch equations. In fact both quantum trajectories and Bloch 
equations are possible and equivalent ways to describe the time evolution of an ensemble 
of fluorescing atoms, but the former is also easy to apply to the emission behavior of a 
single atom. 

As outlined above the photon measurements are taken at rapidly repeated discrete 
times. On a coarse-grained time-scale the paths can be regarded as continuous, and the 
quantum jump approach can be considered as a simple model for continuous measure- 
ments. The photon operators S and J for the photon counting distributions of Sections 3 
in general not only satisfy all the requirements of the axiomatic continuous-measurement 
theory of Ref. |TD[, but are explicitly given, in contrast to the axiomatic theory where 
they have to be more or less guessed. In some cases positivity problems for S and J can 
arise, as discussed further below. 

In Ref. [g \\PoU{t,Q)\%h)\ijW, the no-photon probability at time t was calculated, 
i.e. = II (Op/i|f/(t, 0)|Op/i)|^/') IP, which is, in principle, different from the probability of 
finding no photon until time t, since in the latter case one has to measure in between. 
Interestingly, though, Ref. ^ finds for the no-photon probability at time t the same 
expressions we do for the probability until time t, except for their use of p ■ A coupling 



2^ . This seems to indicate that the reductions used in Eq. (|^) are closely related to the 



Markov assumption and can be replaced by it. Ref. derives photon counting statistics 
by an approach which can be regarded as an alternative to that of Ref. 0, ^ and which 
uses the projector formalism together with the Markov property, considering only finite 
macroscopic time intervals. Interruptions by numerous hypothetical measurements are 
not needed, and neither is the notion of continuous measurements. 

The approximations used in the previous sections are the same as those employed in 



the usual derivations of the Bloch equations |]I8|, P], in particular the use of the Markov 
property, i.e. the fact that the correlation function fi;(r) is very narrowly peaked and drops 
off rapidly. This is a standard assumption in the derivation of Bloch equations [l^, W . 
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The integration of the approximated time derivatives over t' from ti to tj+i in Eq. ([T3|) 
introduces an error for small t' — ti, just as in the derivation Bloch equations and 
a similar error occurs in the reset operator. The omission of the principal values in Eqs. 
(P!^) and ( pUP is not necessary. If one retains them, as for example in Ref. |TB|, only 
the form of the Tiji^s will be changed. The conditional Hamiltonian will then acquire 
an additional term through F, and Eq. ([T8|) should be written in terms of the hermitian 
(real) part Tj- and antihermitian (imaginary) part Fi as 

i^cond = HA{t) + Fi - zFr . (40) 

The imaginary part of F then contributes to the level shifts, but the effect of its possible 
non-diagonal parts deserves closer examination, both for the quantum jump approach and 
for the Bloch equations. 



As seen from Eqs. ( p5[) and (|^), the reset superoperator must preserve positivity, 
i.e. must map positive operators onto positive operators. The question is whether, after 
the approximations used in its derivation, this property still holds. Without explicit 
calculation of the principal value terms in Tijirn little can be said. However, if one drops 
the principal value terms, as in the Bloch equations of Ref. [H, then positivity may be 



lost. In fact, for parallel transition dipole moments the reset matrix of the A system 
calculated from Eq. ( pOf ) then has a tiny negative part for small level separation. 

Similarly, if the transition to a continuous coarse-grained time scale for the conditional 
no-photon time development is possible, then Eq. (^) for the probability density of the 
first photon shows that the hermitian (real) part F^ of F should be a positive operator. 
Without explicit calculation of the principle value part little can be said, but if these are 
omitted then positivity of F^ is lost in some cases when F has non-diagonal elements, e.g. 
for a V system with parallel transition dipole moments. 

These deviations from positivity are probably small in practical applications. But it 
would be preferable to work with expressions which do respect positivity. One way to 
arrive at such expressions is to note that in Eq. ( |30D for J one has, with omission of the 
principal value terms in Tijim, 



D,,-D,^i(4 + ^L) • (41) 



Now, if one replaces the arithmetic mean ^ iuof- + cuf^j by the geometric mean xiuji^uj-i 



'ij ^^Im) '^y "^^^ gcvjiiictin. mean y ujfjUf^ 

then it is easily seen that positivity of the resulting reset operator is automatic. 



Since ti J Pa gives the probability density for the first photon, one has, by Eq. 

tr( Jp^) = tr(F + V*)pA (42) 

if the coarse-grained time scale can be used. Once one has obtained a positivity preserving 
J one can use Eq. (|4^ ) to find a positive V^. A systematic procedure to obtain a positivity 
preserving expression for the reset operator would be desirable. 

In view of the close connection of the quantum jump approach with the Bloch equations 
it stands to reason that the possible non-positivity of the reset matrix and the damping 
operator could have its counterpart in the Bloch equations. Indeed, at least if one omits 
the principal value parts, as in Ref. this can happen for very short times. E.g., for 
the A system with no external driving field and with parallel transition dipole moments 
one can start initially in the upper state which then, under the time evolution under 
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the Bloch equations, develops a small negative part for times of the order of the inverse 
optical frequencies or less, but regains positivity for times larger than that. We have also 
found that a similar phenomenon happens for a V system with parallel transition dipole 
moments. Positivity in the Bloch equations is ensured if one make s the replacement of 
the arithmetic mean in Eq. (^Tj) by the geometric mean. Although this may mean in 
some cases a substantial change of some coefficients in the Bloch equations, the effect on 
the behavior of the solutions is expected to be minimal. A similar non-positivity has been 
found in Ref. 

This seems to be similar to the effect of the rotating-wave approximation. Keeping 
or not keeping the rapidly rotating terms makes some coefficients considerably different, 
but has little effect on the solutions. For the same reason rapidly oscillatory terms arising 
through contributions from Tijj£ for large Uji can be omitted. We have kept them here 
only to exhibit the symmetry of the formulas. 

The derivations in this paper were based on the E ■ D coupling. If one uses the p • A 
coupling the results are the same except that the expression for the Tijim are slightly 
different. If one again replaces the arithmetic by the geometric mean the results become 
the same. 

In conclusion one may say that the degree of reliability of our results for the conditional 
(reduced) Hamiltonian and the reset operator is the same as that of the corresponding 
Bloch equations. Because the times involves are very short - in fact of the order of the 
correlation time used in the Markov approximation - and the non-positivity very small, 
there are probably no practical consequences. 
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